function Uout = hamil(Uin,np,fieldeps,fieldom,jx,jz,E0)

%        Hamiltonian operatioon 

%        phi = H psi

        Uout = 0;
        co= zeros(2);co(1,2)=abs(fieldeps);co(2,1)=abs(fieldeps');
        Uout = (co+diag(E0))*Uin; 
        end